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The nature of the effective spin Hamiltonian and magnetic order in the honeycomb iridates is 
explored by considering a trigonal crystal field effect and spin-orbit coupling. Starting from a 
Hubbard model, an effective spin Hamiltonian is derived in terms of an emergent pseudo-spin- 1/2 
moment in the limit of large trigonal distortions and spin-orbit coupling. The present pseudo-spins 
arise from a spin-orbital locking and are different from the j e g = 1/2 moments that are obtained 
when the spin-orbit coupling dominates and trigonal distortions are neglected. The resulting spin 
Hamiltonian is anisotropic and frustrated by further neighbour interactions. Mean field theory 
suggests a ground state with 4-sublattice zig-zag magnetic order in a parameter regime that can 
be relevant to the honeycomb iridate compound Na2lrC>3, where similar magnetic ground state has 
recently been observed. Various properties of the phase, the spin-wave spectrum and experimental 
consequences are discussed. The present approach contrasts with the recent proposals to understand 
iridate compounds starting from the strong spin-orbit coupling limit and neglecting non-cubic lattice 
distortions. 



I. INTRODUCTION 



Interplay between strong spin-orbit (SO) coupling and electron-electron interaction in correlated electron systems 
has been a recent subject of intensive study 1-24 . In particular, 5d transition metal (e.g. Iridium (Ir) or Osmium (Os)) 
oxides are regarded as ideal playgrounds for observing such cooperative effects 1-24 . Compared to 3d transition metal 
oxides, the repulsive Coulomb energy scale in these systems is reduced by the much larger extent of 5d orbitals, while 
the SO coupling is enhanced due to high atomic number (Z = 77 for Ir and z = 76 for Os). Moreover, owing to the 
extended hd orbitals, these systems are very sensitive to the crystal fields. As a result, the energy scales mentioned 
above often become comparable to each other, leading to a variety of competing phases. Precisely for this reason, 
one expects to see newer emergent quantum phases in such systems. Indeed, there have been several theoretical 
proposals, in context of concrete experimental examples, for spin liquids 1 ' 11,12 ' 14,16 ' 18,19 , topological insulators 10,19 ' 20 , 
Wcyl semimetals 21,22 , novel magnetically ordered Mott insulators 2-6 ' 23,24 and other related phases 8,9 ' 23 in Iridium and 
Osmium oxides. 

A typical situation in the iridates consist of Ir +4 atoms sitting in the octahedral crystal field of a chalcogcn, typically 
oxygen or sulphur 1-4,7 . This octahedral crystal field splits the five hd orbitals of Ir into the doubly degenerate e g 
orbitals and the triply degenerate t^g orbitals (each orbital has a further two- fold spin degeneracy). The e g orbitals are 
higher in energy with the energy difference being approximately 3 eV. There are 5 electrons in the outermost 5d shell 
of Ir +4 which occupy the low lying t-2 g orbitals and the low energy physics is effectively described by projecting out 
the empty e g orbitals 25 . A characteristic feature of most of the approaches used to understand these compounds is to 
treat the SO coupling as the strongest interaction at the atomic level; i.e., by considering the effect of extremely strong 
SO coupling for electrons occupying the t2 g orbitals. This decides the nature of the participating atomic orbitals in 
the low energy effective theory In this limit, the orbital angular momentum, projected to the ti g manifold, carries an 
effective orbital angular momentum i e g = 1 with a negative SO coupling constant 2,3,15 . The projected SO coupling 
splits the t2 g manifold into the lower j e g = 3/2 quadruplet and the upper j e g = 1/2 doublet. Out of the five valence 
electrons, four fill up the quadruplet sector leaving the doublet sector half filled. Thus, in the limit of very strong SO 
coupling, the half filled doublet sector emerge as the correct low energy degrees of freedom. Considering the effect 
of coulomb repulsion within a Hubbard model description and performing strong coupling expansion, various spin 
Hamiltonians for j g = 1/2 are then derived within a strong-coupling perturbation theory 11,20 ' 23 . 

In this paper, we, however, consider a different limit where the oxygen octahedra surrounding the Ir +4 ions are 
highly distorted. While the above scenario of half filled j e s = 1/2 orbitals is applicable to undistorted case, as we 
shall see, it breaks down in presence of strong distortions of the octahedra. In particular, we consider the effect of 
trigonal distortions, which may be relevant for some of the iridate systems including the much debated honeycomb 
lattice iridate, Na2lr03. We show that, in this limit, a different "doublet" of orbitals emerge as the low energy degree 
of freedom. This doublet forms a pseudo-spin-1/2 that results from a kind of (physical) spin-orbital locking, so that the 
spin and orbital fluctuations are not separable (as discussed below). We emphasize that this pseudo-spin is different 
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from the j e s = 1/2 doublet discussed above. The spin Hamiltonian for these pseudo-spins (Eq. 6), on a honeycomb 
lattice, admits a 4-sub-lattice zig-zag (fig. I) pattern in a relevant parameter regime. Such magnetic order has been 
recently observed in the experiments 5,28 on Na2lrC>3 and hence our theory may be applicable to this material. 

The distortion of the octahedron surrounding the Ir +4 generates a new energy scale associated with the change in 
the crystal field, which, as we shall see, competes with the SO coupling. Several kinds of distortion may occur, of 
which we consider the trigonal distortions of the octahedron where it is stretched/compressed along the body diagonal 
of the enclosing cube 25 . In the absence SO coupling, such trigonal crystal field splits the t 2g manifold into e' g (with 
two degenerate orbitals e' lg and e' 2g ) and non-degenerate a\ g (again there is an added two-fold spin degeneracy for 
each of these orbitals). The e' g and a\ g levels are respectively occupied by three and two electrons in Ir 4+ . For large 
trigonal distortions, the splitting between them is big and the a\ g orbitals can be projected out. Now, if one adds 
SO coupling, the low energy degrees of freedom is described by a subspace of the e' g orbitals which form an emergent 
pseudo-spin-1/2 doublet out of |e' lg ,|) and \e' 2g |) states, where f and J, represent the physical spin s z — 1/2, —1/2 
(the spins arc quantized along the axis of trigonal distortion). These pseudo-spin-1/2 is different from the j e g = 1/2 
and j e s = 3/2 multiplets in the strong SO coupling limit as discussed above. Notice the (physical)spin-orbital locking 
for the pseudo-spins, as alluded above. 

The two approaches, described in the last two paragraphs, of arriving at the low energy manifold are mutually 
incompatible. This can be seen as follows: In presence of sizeable trigonal distortions the j e ff = 1/2 and j ff = 3/2 
multiplets mix with each other and can no longer serve as good low energy atomic orbitals. This dichotomy becomes 
quite evident in the recent studies of Na2lr03, where, the Jr +4 form a honeycomb lattice. Taking into account the 
strong SO coupling of Ir 4+ in Na2lr03, proposed are a model for a topological insulator in the weak interaction limit 10 
and a Hciscnbcrg-Kitacv (HK) model for a possible spin liquid phase in the strong coupling limit 11 . These proposals 
prompted several experimental 4-6 and theoretical efforts 12,26 to understand the nature of the ground state in this 
material. Subsequently, it was found that Na2lr03 orders magnetically at low temperatures 4 . However, the magnetic 
moments form a "zig-zag" pattern (fig. I) which is not consistent with the ones that would be obtained by adding a 
weak interaction in a topological insulator (canted antiferromagnet) 10 or from a nearest neighbour HK model (spin 
liquid or the so-called stripe antiferromagnet) 11 . While recent studies 27 show that a 'zig-zag' order may be stabilized 
within the HK model by including substantial second and third neighbour antifcrromagnetic interactions, it is hard 
to justify such large further neighbour exchanges without lattice distortions. (An alternate explanation that we do 
not pursue here is significant charge fluctuations which would mean that the compound is close to metal-insulator 
transition. The resistivity data seems to support the fact that this compound is a good insulator. 4 ) If lattice distortions 
are responsible for the significant further neighbour exchanges, then, there would be sizeable distortion of the oxygen 
octahedra, which in turn may invalidate the above j c s = 1/2 picture and thus the basic paradigm of the HK model, 
by mixing the j e g=l/2 and j e ff = 3/2 subspaccs. Recent finite temperature numerical calculations on the HK model 26 
also suggest possible inconsistencies with experiments on Na2lr03. This necessitates the need for a different starting 
point, to explain the magnetic properties of Na2lr03. 

The rest of this paper is organized as follows. In Section II, we derive the effective spin Hamiltonian in limit 
of the large trigonal distortion and large SO coupling. This is done by taking the energy scale associated with 
trigonal distortion to infinity first, followed by that of the SO energy scale. This order of taking the limit gives a 
spin Hamiltonian in terms of emergent pseudo-spin— 1/2, which is different from the HK model. This Hamiltonian 
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is both anisotropic and frustrated. It also has further neighbour interactions, the effect of which are enhanced due 
to anisotropy that makes some of the nearest neighbour bonds weaker. The origin of this anisotropy is trigonal 
distortion. We argue that this limit may be more applicable for the compound Na2lrC>3. Having derived the spin 
Hamiltonian, we calculate the phase diagram and the spin wave spectrum within mean field theory in Sec. III. We see 
that the 'zig-zag' phase occurs in a relevant parameter regime. We also point out the experimental implications of our 
calculations in context of Na2lrC>3. Finally we summarize the results in Sec. IV. The details of various calculations 
are given in various appendices. 



II. THE EFFECTIVE HAMILTONIAN 



In the cubic environment the t 2g orbitals are degenerate when there is no SO coupling. Trigonal distortion due to 
compression or expansion along one of the four C3 axes of IrC>6 octahedra lifts this degeneracy. Although it is possible 
that the axes of trigonal distortions are different in different octahedra 20 , we find that uniform distortions are more 
consistent with the experiments (see below) on Na2lrC>3. Hence we consider uniform trigonal distortion. 



A. The Trigonal Hamiltonian 



Let us denote the axis of this uniform trigonal distortion of the octahedron by the unit vector n = -i=[m, n 2 , 113], 
where n a = ±1. Since there are 2 directions to each of the 4 trigonal axes we may choose a "gauge" to specify n. This 
is done by taking n\n 2 nj, = +1. The Hamiltonian for trigonal distortion, when projected in the t 2g sector, gives 
our chosen gauge) 



2(1 



Hi 



n 3 n 2 
n 3 ni 
ni Tii 



where ^\ = [d yz ,d\ x ,d xy ] and A tr ; is the energy scale for trigonal distortion. The eigenstates are 
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The trigonal distortion splits the t 2g sector into the doubly degenerate e' g and the non-degenerate a\ g with energies 
Atri/3 and — 2A tr i/3 respectively. 

A description based on Hubbard model for the e' g orbitals may be systematically derived starting from the t 2g 
orbitals. This is done in A. This has the following general form 



H' = H 



SO 



e iMa e iM'a' e iM'<j' S-iMcr , 
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where e iM(T is the electron creation operator in the e' g orbital (M = 1, 2) with spin cr(=f, J,); Um-jm' are the effective 
hopping amplitudes within the subspace and U is the effective onsite coulomb's repulsion. We note that the Hund's 
coupling (which arises from the orbital dependence of the Coulomb repulsion) for the t 2g orbitals only renormalizes 
U in this restricted subspace. (See A for details). 



B. The Projected SO coupling 



The SO coupling, when projected in the e' g subspace, yields a block diagonal form (see B for details): 

Hs° = -\n-s iT ?, (4) 
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FIG. 2: The e' states split by the SO coupling. 



where Si is the spin operator at the site i, A « 500 raeF is the SO coupling parameter and = +1(— 1) refers to the 
e i ff ( e 2g) orbital. 

Thus the projected SO interaction acts as a "Zeeman coupling" where the direction of the "magnetic field" is along 
the trigonal axis or opposite to it 13 . Thus it is natural to choose the direction of spin quantization along the axis 
of trigonal distortion. This then gives the active atomic orbitals after incorporating the SO coupling. These active 
orbitals are the Krammer's doublet |e' lg ,4-) and |e2„,t) as shown in Fig. 2. 



C. The Spin Hamiltonian 



Hence the low energy physics may be described by considering only the above atomic orbitals. The starting point 
for the calculations is projection of the Hubbard model (Eq. 3) in the space spanned by the Krammcrs's doublet 
\e' lg ,l) and |e 2 „,t}- The bandwidth of this projected model is narrow and the effect of the Hubbard repulsion is 
important. Indeed it can easily render the system insulating. To capture the magnetic order in this Mott insulator, 
we do a strong-coupling expansion in t/U to get an effective "pseudo-spin" model in terms of the pseudo-spin- 1/2 
operators, 

S a = \e\ P a ab e bl (5) 

where, p a (a = x,y,z) arc the Pauli matrices and a,b — (e' g i', I), (e' g2 ',t)- The "pseudo-spin" Hamiltonian has the 
following form up to the quadratic order: 



H - E<ij) J iP&i ■ + J^Si ■ Sj + £(((«>» J ifSi ■ Sj 

i t( z 1) qz qz i t( z2 ) q z qz i t( z 3) q z q z 

Here (ij), ((ij)) and (((ij)}} refer to summation over first, second and third nearest neighbours (NNs) respectively. 
The different exchange couplings arc given in terms of the underlying parameters of the Hubbard model as 
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where a = 1,2,3 denotes that ij arc first, second or third NNs, respectively and the last expression defines Jij^ ■ 

T^°^ and T>^ a ' are given in terms of the hopping amplitudes (e.g. t^' yz from the overlap of xy and yz orbitals) of 
the t 2g orbitals as (details are given in C). 



(Or 



{{j.VZ\vz . ,xz\xz 

3 .('" + V 

^ \( ( ,xz\xy | ,xy 



,xy;xy 



)] 



+n 2 a:^- 



+ n 3 (i 



yz;xz 



(8) 



TAza) 1 r { ,xz\xy ,xy\xz\ . f.xy.yz ,yz:xy\ , { ,yz\xz ,xz;yz\~\ 

h = irm C ni (*« " *« ) + na (*« " ) • »•> (*« - V )J 



(9) 



FIG. 3: Section of a honeycomb lattice (shaded in yellow). Ir sites (black) are connected by bonds (orange). The green arrow 
is the the [— 1, — 1, 1] direction of trigonal distortion that makes an angle of about 19° with the plane of the lattice pointing 
inside the plane. 

Before moving on to the details of the spin Hamiltonian, we note that, on projecting to the subspace of |e' lff ,4.) and 
| it)) the spin and orbitals are no longer independent. Instead at every site there is a pseudo-spin- 1/2 degree of 
freedom where the spin is locked to the orbital wave function. This, we refer to as spin-orbital locking. 

III. APPLICATION TO Na 2 Ir0 3 

We now apply the above results to the case of Na2lrC>3. The early X-Ray diffraction experiments 4 suggested a a 
monoclinic C2/C structure for the compound and distorted Ir06 octahedra. However, more recent experiments see a 
better match for the X-Ray diffraction data with the space group C2/m. 28,30 They also unambiguously confirm the 
presence of uniform trigonal distortion of the IrOg octahedra. However, the magnitude of such distortion is not clear 
at present. Further, experimental measurements suggest: (1) The magnetic transition occurs at Tn = 15K while the 
Curie- Weiss temperature is about &cw ~ — 116if . This indicates presence of frustration. (2) The high temperature 
magnetic susceptibility is anisotropic; the in-plane and out-of-plane susceptibilities are different. This may be due to 
a trigonal distortion of the 1x06 octahedra 4 . (3) The magnetic specific heat is suppressed at low temperatures. 4 (4) 
Recent resonant X-ray scattering experiment 5 suggests that the magnetic order is collinear and have a 4-site unit cell. 
(5) The magnetic moments have a large projection on the a-axis of the monoclinic crystal 5 . (6) A combination of 
these experimental findings and density functional theory (DFT) calculations strongly suggest that a 'zig-zag' pattern 
for the magnetic moments, as shown in Fig. I in the ground state 5 , which has since been verified independently by 
two groups using Neutron scattering 28,30 . 

Taking these phcnomenological suggestions, wc try to apply the above calculations to the case of Na2lr03. At 
the outset, we must note that, in the above derivation of the spin Hamiltonian we have assumed that the trigonal 
distortion to be the largest energy scale followed by the SO coupling. While this extreme limit of projecting out 
the a\ g orbitals most likely is not true for Na2lr03. However, we expect the real ground state to be adiabatically 
connected to this limit. With this in mind, we now consider the case of Na 2 Ir03. 

Clearly, the exchanges (Eq. 7) depend both on the direction of the bond and the direction of the trigonal distortion. 
So it is important to ask about the direction of the latter. Comparing the crystallographic axes of Na2lr03, we find 
that the direction [1,1,1] is perpendicular to the honeycomb plane while the other three directions make an acute 
angle to it. In the monoclinic C2/m structure, uniform trigonal distortion in these four directions may not cost the 
same energy. In experiments 5 , the moments are seen to point along the a-axis of the monoclinic crystal which is 
parallel to the honeycomb plane. This, along with the fact that the magnetic moment in our model is in the direction 
of n (explained below) seems to suggest that n = ^=[— 1, — 1, 1] is chosen in the compound (see Fig. 3). In the absence 
of a better theoretical understanding of the direction of the trigonal distortion, we take this as an input from the 
experiments. 

To identity different hopping paths (both direct and indirect), we consider various overlaps (see C) and find, while 
j(3z) = o ; j( lz ) ^ j( 2z ) ^ are approximately (spatially-)isotropic and antiferromagnetic. For the exchanges of the 
Hciscnbcrg terms, both and are antiferromagnetic and isotropic (both of them result from indirect hopping 
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FIG. 4: Mean field phase diagram for Eq. 6. The two axes are: xo = y^', Ho = ■jpy> where J^(J^) are related to the 
strong(weak) NN exchange and j' 2 - 1 is the 2 nd and 3 rd neighbour exchange (see Eq. 7). We take the Ising anisotropy to be 5% 
of J* 1 '. Note that, due to Ising anisotropies, one has zig-zag order at j/o = 0. 

mediated by the Na s-orbitals and are expected to be comparable). For the NN Heisenberg exchanges, the couplings 
are antifcrromagnetic, but, much more spatially anisotropic. We find that for the chosen direction of the trigonal 
distortion, the coupling along one of the NN exchanges (J*- 1 -*) (viz b\ in Fig. 3) is different from the other two 
neighbours (J^)(&2 and 63 in Fig. 3). 

A. Mean-Field Theory and Magnetic Order 

We now consider the mean field phase diagram for the above anisotropic spin Hamiltonian. For being the 
largest energy scale, the classical ground state for the model can be calculated within mean-field theory as a function 
of xq = JW j 'JW and yo = J^/JW (we have taken — J^). A representative mean-field phase diagram is shown 
in Fig. 4. It shows a region of the parameter-space where the zig-zag order is stabilized 29 . The effect of the Ising 
anisotropies j' lz ) and j( 2z ) is to pin the magnetic ordering along the z-direction of the pseudo-spin quantization which 
is also the direction of the trigonal distortion n. They also gap out any Goldstone mode that arises from the ordering 
of the pseudo-spins. The latter results in the exponential suppression of the specific heat at low temperatures. The 
other competing phase with a collinear order is the regular two-sublattice Neel phase. 

The nature of the ground states may be understood from the following arguments. In the presence of the n in 
[—1,-1,1] direction, the NN exchange coupling becomes anisotropic. When it is strong in one direction (J^ 1 * 1 ) and 
weak in two other directions (J*- 1 -*), for the bonds where the NN coupling becomes weak, the effects of the small 
second and third neighbour interactions become significant. Since the latter interactions are antiferromagnetic, they 
prefer anti-parallel alignment of the spins. As there are more second and third neighbours, their cumulative effect can 
be much stronger. This naturally leads to the zig-zag state. The NN antifcrromagnetic interactions on the weaker 
bonds compete with the antifcrromagnetic second and third neighbour interactions and frustrates the magnet. This 
suppresses the magnetic ordering temperature far below the Curie- Weiss temperature. 

B. The spectrum for Spin-orbital waves 

The low energy excitations about this magnetically ordered zig-zag state are gapped spin-orbital waves. Signatures 
of such excitations may be seen in future resonant X-Ray scattering experiments. It is important to note that this 
"pseudo-spin" waves actually contain both orbital and the spin components due to the spin-orbital locking. 

We calculate the dispersion of such spin-orbital waves to quadratic order using the well-known Holstcin-Primakoff 
methods. The details are discussed in D. A representative spin wave spectrum in the zig-zag phase is shown in Fig. 5(a) 
and 5(b). The spectrum is gapped and the bottom of the spin- wave dispersion has some characteristic momentum 
dependence. 

C. Experimental Implications 

Apart from the already discussed exponential suppression of low temperature magnetic specific heat, the above 
calculation predicts an interesting feature in the magnetic susceptibility The relation between the magnetic moment 
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FIG. 5: The "pseudo-spin" wave spectrum (contours of both the bands are shown in (a) and a section is shown in (b)). The 
values used for the parameters are same as that used for the calculation of the mean field phase diagram (Fig. 4). We note 
that, as expected, the spectrum is gapped. 

and the pseudo-spins is 

Mi = -4»BnS?, (10) 

where [Ib is the Bohr magneton. This follows from the twin facts that, in e' g subspacc, the angular momentum 
transverse to n is quenched and the spins are locked to the orbitals with the axis of quantization being n in our 
pseudo-spin sector (see E). Thus, the magnetization is sensitive to the z-component of the pseudo-spin (the direction 
of which is shown in Fig. I). Indeed the magnetization has the largest projection along the a- axis of the monoclinic 
crystal. This was seen in experiments 5 and was the motivation for choosing the [—1,-1,1] direction for the trigonal 
distortion. Along two other axes [—1, 1, —1] and [1, —1, —1], a large component of in-plane magnetization exists, but in 
different directions. Finally the direction [1,1,1] is perpendicular to the honeycomb plane and leads to magnetization 
in the same [1,1,1] direction. While this does not appear to be the case for Na2lr03, this may be more relevant for the 
less-distorted compound Li2lr03 (see below). Eq. 10 suggests that the magnetic susceptibility is highly anisotropic 
and depends on the cosine of the angle between the direction of magnetic field and n. Indeed signatures of such 
anisotropy have been already seen in experiments 4 . We emphasize that within this picture, the in- plane susceptibility 
also varies with the direction of the magnetic field. So the ratio X±/x\\ can be lesser or greater than 1. The current 
experiments 4 does not tell the in-plane direction of magnetic field and hence we cannot comment on the ratio presently. 
However, the above picture is strictly based on atomic orbitals. One generally expects that there is also hybridization 
of the Ir d-orbitals with the oxygen p-orbitals. Such hybridization will contribute to a non-zero isotropic component 
to the susceptibility 24 . Also, as remarked earlier, in the actual compound, the SO coupling scaling may not be very 
small compared to the trigonal distortion limit scale. Additional perturbation coming from the mixing with the a\ g 
orbitals will also contribute to decrease the anisotropy of the susceptibility. 
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IV. SUMMARY AND CONCLUSION 



In this paper, we have studied the effect of trigonal distortion and SO coupling and applied it to the case of the 
honeycomb lattice compound Na2ir03. Wc find that, in the limit of large trigonal distortion and SO coupling, a 
pseudo-spin- 1/2 degree of freedom emerges. Low energy Hamiltonian, in terms of this pseudo-spin gives a 'zig-zag' 
magnetic order as seen in the recent experiments on Na2lr03. We have also calculated the low energy spin- wave 
spectrum and elucidated various properties of the compound that has been observed in experiments. The pseudo-spin 
couples the physical spin and the orbitals in a non-trivial manner, signatures of this may be seen in future inelastic 
X-ray resonance experiments probing the low energy excitations. 

While very recent experiments 28,30 clearly indicate presence of trigonal distortions, their magnitude is yet not 
confirmed. On the other hand, the only available numerical estimate of the energy scale for trigonal distortion comes 
from the DFT calculations by Jin et al. 13 (based on C2/C structure). It suggests A tr i ~ 600 raeV . While, it is not 
clear if such a large value is in confomity with the experiments, at present, the detection of trigonal distortion in 
experiments is highly encouraging from the perspective of the present calculations. 

In these lights of the above calculations, it is tempting to predict the case of Li2lr03 where recent experiments 
suggest a more isotropic honeycomb lattice 6,31 . A possibility is that sizeable trigonal distortion is also present in 
Li2lr03 (so that the above discussion holds), but, the axis is perpendicular to the plane. What may be the fallouts in 
such a case ? Our present analysis would then suggest that the antiferromagnetic exchanges are isotropic and equally 
strong for the three NNs. This would develop 2-sublatticc Necl order in the pseudo-spins with the magnetic moments 
being perpendicular to the plane. Also the further neighbour exchanges are rather weak (compared to Na2ir03) and 
hence frustration is quite small. Indeed recent experiments see ordering very close to the Curie- Weiss temperature, 
the later being calculated from the high temperature magnetic susceptibility data 6,31 . However, present experiments 
do not rule out the possibility of small or no trigonal distortions in Li2lr03, in which case the limit of HK model 11 ' 26 
may be appropriate. 
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Appendix A: The microscopic model for Na2lr03 



The generic Hubbard model (for the t 2g orbitals) including the trigonal distortions, Hund's coupling and the SO 
coupling is 



+ 2 Ei Emm' Saa' U m m' d ima d im , a , dr ml t a i di ma ■ 



(Al) 



Here m, m' — yz, xz, xy and a =t, 4- an d H t ^ is given by Eq. 1. We note that the hopping is diagonal in spin space 
and in the cubic harmonic basis all hopping arc real. Also, the hopping contain both the direct and indirect (through 
Oxygen and Sodium) paths. We have taken Hund's coupling into account through U mm i, though this is expected to 
be small in 5d transition metals. To a very good approximation the form of U mm i is given by 



Un 



U U () -Jh Uo-Jh 
Uo — Jh Uq Uo — Jh 
Uo-Jh U -J h U 



(A2) 



where the basis is given, as before, by VfJ = [d^ z , d\ x , d\. y \. Uq and Jh are the intra orbital Coulomb repulsion and 
Hund's coupling term respectively. 



The transformation between the operators in the trigonal basis, & 
[dy z: d\ x , d xy \ , is given by ^ m = T„ 1: m&m- The transformation matrix is given by 

T m ,M 



t 't 't 

l 9 ' e l ff ' e 2 ff 



and t 2g basis, & 
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(A3) 
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The transformations for the hopping amplitudes and repulsion term are then given by 



HM;jM> — 2sm,m' 1 m,Mij 1 m',M'i 

U Ml hh = J2 m ,m' U mm' (TmM 1 T mM 1 ) ( T m'M 2 T m'M 2 ) ■ (A4) 



Notice that there are contributions to i™'" 1 from both direct and indirect exchanges for the first, second and third 
neighbours, as confirmed from the DFT calculations by H. Jin et al. 13 . These show that there are contributions from 
both direct and indirect hoppings for the first, second and third nearest neighbours. Projecting them into the e' g 
orbitals we get the effective hopping amplitudes which are then used in Eq. 3. As for the Coulomb repulsion term, 
we find that it has the following form 



Um x m 2 = U 



1 1 1 
1 1 1 
1 1 1 



(A5) 



where U = Uq — 2J#/3. This form is then used in Eq. 3. The reason for this special form of Um 1 m 2 ncs m the fact 
that the e' g orbitals have equal weight of the three t 2g orbitals (see the wave functions in Eq. 2). 



Appendix B: Projection of Spin-Orbit coupling to the e' g subspace 



The SO coupling, when projected to the t- 2g orbitals give 

Hs6 = -Af- s, (Bl) 

where I is a I = 1 angular momentum operator. We can re-write the t 2g cubic harmonics in terms of the spherical 
harmonics of the effective I = 1 angular momentum operator. These are given by: 

&,) =^[|1,-1)-|1,+1)]; 
\d zx ) =^[|1,-1) + |1,-1)]; 

\d xy ) =|1,0) (B2) 

The projector for the e' g space is: P e s = |ei)(e^| + | e^) (e^ | - It turns out that I - sis block diagonal in this subspace. 
Hence 

f- s= KXellf- s\e\){e\\ + |e' 2 )(e' 2 |f- S >' 2 )(e' 2 | (B3) 
Making the "gauge" choice we get 

(e[\l ■ s\e[) = n ■ s; (e' 2 \l ■ s\e' 2 ) = -h ■ s (B4) 



Appendix C: The hopping parameters 



1. Nearest neighbours 



The nearest neighbours are shown in Fig. 6(a). There are two different processes contributing to the hopping.: 1) 
the direct hopping between the Ir atoms and 2) the indirect hopping between the Ir atoms mediated by the oxygen 
atoms. In presence of the trigonal distortion which has a component along the honeycomb plane (like in this case 
[— 1, — 1, 1]) the magnitudes of the different hopping parameters arc different in different directions (for both direct 
and indirect hopping). The results are shown in Table I. 

We shall make an approximation here. We shall leave out the directional dependence of the magnitudes on the 
direction. The argument is that the essential directional dependence due to the trigonal distortion has been taken 
care of by the parameter Ai. When the DFT 13 results are used to find the tight-binding parameters 29 , it is found 
that (they use Ai = 0) (here t^di and t^i are direct hopping and to is the indirect hopping respectively.) tddi = 
-0.5 eV;t dd2 = 0.15 eV;t Q = 0.25 eV. 
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FIG. 6: The 3 nearest neighbours (a), six 2 nd nearest neighbours (b) and three 3 rd nearest neighbours (c) of the central site. 
The nomenclature has been used to label the hoppings. 



(a)NN:t om . 6l , 



\ m 
m' \ 


dxy 


dyz 


d zx 


dxy 


tddi{bi) 










tdd2(bi) 


-t dd2 {bi) +to(6i) + Ai(6i) 


dzx 




-tdda(bi) + to(bi) - Ai(bi) 


tdds{bi) 



(b)NN:t am]b2r 



\ m 
m 1 \ 


dxy 


dyz 


dzx 


dxy 


tdd2{b 2 ) 




-t d d2(b 2 ) + t (b 2 ) + Ai(b 2 ) 


dyz 








d zx 


-tdda(6a) +*o(62) - Ai(6 2 ) 




tdd2{b 2 ) 



(c)NN: t am . b3T 



\ m 
m' \ 


dxy 


dyz 


dzx 


dxy 


^2(63) 


-tdd2(b :i ) +t (b 3 ) + Ai(6 3 ) 




dyz 


-t d d2(b 3 ) +t (b 3 ) - Ai(6 3 ) 






d Z x 






tddi{bz) 



TABLE I: The hopping paths (both direct and indirect) in the t 2g basis. 



Performing the transformation to the e' basis, we have 



and 



Hence, 



T 



(02) _ 1 



ab-> 



T, 



(03) 



ah,' 



{tdd2 - to)n 3 ] , 
{tdd2 - to)ni] , 

I [tddl + 2tdd2 + (*dd2 ~ to)n 2 ] ■ 



— 3 [tddl + 2tdd2 



T (zl) _ _ Ai 



T 



(zl) _ 



T, 



(zl) 



V3 



=n 2 



t(0) 



J, 



4 

3! ' 



a b 

(0) _ , 



3fr 



i [tddi + 2t dd2 + (i drf2 - i )n 3 ] 2 - (A1) 2 
3 [irfdi + 2t dd2 + (idrf 2 - to)m] 2 - (A1) 2 
I [*ddi + 2t dd2 + (t dd2 - t )n 2 ] 2 - (Ax) 2 



(CI) 



(C2) 



(C3) 
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(12) 

ah i 
(U) 
ab 2 
(I2) 



8(Ai)- 
3(7 

8(Ai) s 



3(7 
8(AQ 



-abi 3(7 

where we have taken the direction of the trigonal distortion is taken to be uniform. 



(C4) 



2. Second nearest neighbour 

These are shown in Fig. 6(b). These indirect hoppings are mediated by the Na atoms. In general, in presence of the 
trigonal distortion in the [—1,-1,1] direction, the magnitude of the hopping amplitudes are also direction dependent. 
However, since the magnitudes themselves are expected to be small we shall neglect such directional dependence in 
the magnitudes. The result is summarized in Table II. 

(a)NNN.f am . aim //t am . a ^ m / (b)NNN.t am . (l2m //t flm . (I ^ m / (c)NNN. iam;a 3 m' /^am;aQm' 
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dy Z 
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dxy 




t 2 + A 2 




dyz 


t 2 - A 2 






dzx 









\ m 
m' \ 


dxy 


dyz 


dzx 


d X y 






t 2 + A 2 


dyz 








dzx 


t-2 - A 2 







m' \ 


dxy 


dyz 


dzx 


dxy 








dyz 






t 2 + A 2 


dzx 




t 2 ~A 2 





TABLE II: Hopping paths for the second nearest neighbours 



So for the e' g basis, we have 



T (02) 
a,ai 


_ T (02) 
a, a<i 


= -|n 2 , 


T (z2) 


_ T (z2) _ 


A 2 

-7§ n2 


T (02) 
a,d2 


_ T (02) 
a,a 5 


= -|m, 


rp(z2) 

± a,a 2 


T (z2) 
— -L a,a 5 — 


A, 

"7t ni 


T (02) 
a,a 3 


- T (02) 


= -|n 3 , 


T (z2) 


_ T (z2) _ 


A, 



For example, tight binding fit of the DFT data uses only t2 and finds £2 ~ —0.075 eV ' . Therefore we have 



a,a a w 



(t2? 



~ (A 2 ) 5 



r(2z) _ 



8(A 2 ) S 
3U 



(C5) 



(C6) 



3. Third nearest neighbour 

The third nearest neighbours are listed in Fig. 6(c). The hopping to the third nearest neighbour is mediated by the 
Na atoms. Again we shall neglect the directional dependence and take these to be in the magnitudes of the hoping 
amplitudes. The result is summarized in table III. 



(a)NNNN:t 6 / 
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*s(6i) 






dyz 








dzx 









(b)NNNN:t am b / 
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Hb' 2 ) 
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(c)NNNN: t am . b . 



\ 711 

m' \ 


dxy 




dzx 


dxy 








dyz 








dzx 






t 3 (b' 3 ) 



TABLE III: The hoppings for the third nearest neighbours 
Tight-binding fit to the DFT results 13,29 indeed show that this hopping energy scale is of the order of 



ts(b' a ) = t n Ps -0.075 eV 



Therefore we have: 



± ab' 



(03) _ in 



1 1 nh' 



(23) 



(C7) 
(C8) 
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or. 



J; 



(3) 



9U 



T (3z) 



(C9) 



Appendix D: Spin Wave Spectrum 



To calculate the spin wave spectrum for the zig-zag state we use the usual Holstein-Primakoff method suited to 
collinear ordering which may alternate in direction. More precisely we introduce: 



for one direction and 



S z = S - a^a; S+ = V2Sa; S~ = V2Sa^ 



S z = -S + a) a; S+ = V^Sa^; S~ = V2Sa 



(Dl) 



(D2) 



for the other direction. Since there are 4 sites per unit cell (refer Fig. 1(a) of the main text) the quadratic Hamiltonian 
is a 8 x 8 matrix given by: 

H Q = H cl + H sp , (D3) 
where H c i is the classical part dealt in the previous section. The spin wave Hamiltonian has the following form 



H sp = |^k# k K k * k + H s 



(D4) 



a lc i' a iL 2' a iL 31 a iL 4' a -k,i 5 a -k.2, <2-k,3, Q-k,4 1 (the subscript 1,2,3,4 refers to the four sites in the unit 



Here w u 

cell as shown in Fig. 1(a) of the main text) and 



H s = -f [(1 - 2x + by) - (26 2 - St)} N c 



ell • 



where N ce ii is the number of unit cells and 



A^ S k 
Bt A u 



(D5) 
(D6) 



Xk o r? k 

Xk 0k o 

o <f>i x k o 

r& xk 



& Pk 

a o o Pk 

pi o o Ck 

L o pi a o 



where, 



Xk = (252- 6 1 ) + (l-2x + 5y + ycosk x ); 



'/k 

0k 
PU 

Now following usual methods we diagonalize 



= xe tky (l + e lfejc ) ; 
= aj(l + e**-) ; 
= (1 + 2ycosk x + ye~ tk y); 
= y{\ + e lk *)(l + e~ lk y). 



-Bt -Ak 



(D7) 



(D8) 
(D9) 
(D10) 
(Dll) 
(D12) 



(D13) 



to get the spin wave spectrum as plotted in Fig. 3(a) and 3(b) of the main text. 
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Appendix E: Projection of Zeeman term in the t 2g and {\e' lg l), \e' 2g t}} subspaces. 

The Zeeman coupling term, when projected to the t2 g space, gives 

= Li B f-f+ 2s) • B. (El) 

Thus the magnetization after projection is given by: 

A?' 2j = fi B (-T+ 2s) (E2) 

This when projected to the subspace |e' 1; f) and \e!%, l) gives (using the Block diagonal property of the orbital angular 
momentum as above): 

H z = 4:fi B S z h ■ B (E3) 

where S (note that this is in upper case compared to the physical spin written lower case) is the emergent pseudo- 
spin-1/2 per site. This is the emergent degree of freedom at low energies. Clearly, the magnetic moment is then given 
by Eq. 7 of the main text. 
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